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Modulus  of  elasticity  of  the  blade  (lb/in2) 
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Complex  modulus  of  elastomer  (Ib/in  ) 
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Real  part  of  complex  modulus  of  elastomer  (lb/in  ) 

Imaginary  part  of  complex  modulus  of  elastomer 
( lb/in" ) 
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Area  moment  of  Inertia  of  the  blade  (in  ) 
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Mass  moment  of  inertia  lb-sec  -in 
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NOTE:  Computer  program  in  Appendix  C  restates  some  of  the 
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compressor  blade  as  a  means  of  damping  vibratory  oscillations. 
The  avenues  to  this  exploration  were  an  experimental  and  an 
analytical  model  of  the  blade  along  with  a  comparison  and 
analysis  of  these  models. 
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Th.e  derivation  of  an  equation  of  motion  satisfying  par¬ 
ticular  boundary  conditions  was  the  main  objective  of  the 
analytical  section.  The  turbi.ne/corapressor  blade  was  modeled 
as  a  modified  cantilever  beam  with  a  torsional  spring  and 
viscous  damper  at  the  root  to  represent  the  damping  material, 
and  a  tic  mass  with  some  added  stiffness  at  the  tic  to 
account  for  the  shaker  stio'fness  and  the  mass  of  the  shaker, 
force  gauge,  linkage  and  accelerometer.  The  solution  of  the 
resulting  eigenvalue  problem  was  used  to  prediot  the  resonant 
frequencies  of  the  blade. 


mi  .n,  '  r  o  'rr  •  c.  c?  ->  —  v->  -5  •  o  —  a  ~~  r-,  -j  tc  ^  o 

—  .  «•  W  ‘  -  -  -*  —  — .  .*  -  .  -1 ,J  b  V  • »  —  O  V  ^  •  lU  U  ^  ^  ^  ^  ^  ^  ~  —  .  .  .  ^  - 

using  a  commercially  available  damping  material.  Each  vis¬ 
coelastic  material  has  certain  frequency  and  temperature 
ranges  where  losses  reach  a  maximum.  Consequently,  a  certain 
material  can  be  chosen  to  fit  a  particular  application.  If 
viscoelastic  damping  at  the  blade  root  proves  feasible, 
further  studies  employing  materials  usable  at  engine  oper¬ 
ating  temperatures  should  be  pursued.  It  is  already  well 
known  that  vitreous  materials  have  high  loss  factors  that 
can  be  employed  to  considerably  reduce  resonant  vibratory 
energy  (Ref  1).  Such,  materials  may  be  desirable  in  this 
application. 
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5  shows  the  spring  beam  that  was  calibrated  to  provide  the 
contact  pressure  that  also  acted  to  simulate  the  centrifugal 
forces  developed  in  a  rotating  jet  engine.  The  calibration 
of  this  beam  is  outlined  in  Appendix  A. 

The  viscoelastic  material  used  in  the  experiment  was 
ISD  112.  It  was  chosen  because  it  has  peak  damping  proper- 

-i 

ties  at  ambient  temperatures  and  at  frequencies  between  ID*1 


~  I 

and  lO"5  cycles  per  second  which  were  the  test  conditions.  i 
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BLADE  CONFIGURATION  AND  DIMENSIONS 
FIGURE  4 


Table  1,  an  extract  from  reference  8,  lists  the  actual 
properties  of  the  material  used  in  the  tests. 

Table  2  delineates  the  properties  of  the  blade  and  the 
dovetail  root. 

The  mass  and  stiffness  of  the  shaker  were  determined 
experimentally  as  delineated  in  Appendix  3. 

Test  Set-Up 

A  schematic  diagram  of  the  test  set-up  is  shown  in  Figure 
1.  The  shaker  exciter  was  a  MB  Vibration  test  equipment  model 
T122133-  This  was  used  to  vary  the  frequency  and  amplitude 
of  the  sinusoidal  signal  from  the  shaker.  The  frequency  was 
monitored  using  a  Computer  Measurements  frequency  counter, 
model  725 C.  The  shaker,  an  M3  model  SC,  was  used  tc  simulate 
a  sinusoidal  tip  loading  that  would  represent  the  periodic 
force  on  a  turbine/compressor  blade.  A  ?C3  force  transducer 
model  2G8AC3  serial  number  1924  was  mounted  in  series  with  the 
vibrator  linkage.  For  this  device,  the  vendor  supplied  the  - 
calibration  factor  of  11.6  millivolts  per  pound.  The  trans¬ 
ducer  output  (see  figures  2  and  u)  was  then  amplified  by  the 
variable  gain  PC3  amplifier  model  48CC06  power  unit.  The 
accelerometer  (see  figures  3  and  4)  was  a  PCB  model  303A  serial 
number  2066  with  a  calibration  factor  of  14.3  millivolts  per 
g  that  was  also  supplied  by  the  vendor.  It  toe  was  amplified  ~ 
by  a  similar  PCB  power  unit.  Both  the  force  gauge  and  the 
accelerometer  outputs  were  measured  using  a  Tektronix  oscillo¬ 
scope  model  535A.  The  copper  constantan  thermocouple  was  spot 
welded  to  the  blade  mount  in  the  vicinity  of  the  elastomer 
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37.3 

284.3 
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E+05 
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1.1351 

37.3 
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BLADE  AND  DOVETAIL  PROPERTIES 


Blade  and  Dovetail  Material 

Ir.conel-X 

Modulus  of  Elasticity 

31(lo6)  ib/in2 

Blade  area  Moment  of  Inertia 

0.001373  in4 

Blade  Mass 

0.003905  lb-sec2-in 

Dovetail  Mass 

0.000719  lb-sec2-in 

Total  Mass 

0.004624  lb-sec2-in 

31ade  Mass  Moment  of  Inertia 

about  point  3  (see  figure’ 4) 

0.1402  lb-sec2-in 

Dovetail  Mass  Moment  of  Inertia 

about  point  0 

0.00001975  lb-sec2-in 

Total  Mass  .Moment  of  Inertia 

about  point  0 

0 . 14021975 lb-sec2-in 

TABLE  2 
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Test  Procedure 

For  each  frequency  sweep,  the  displacement  at  the  tip  of 
the  blade  was  held  constant  by  setting  the  appropriate  ac¬ 
celeration  as  determined  in  step  5.  The  strain  gauges  on  the 
spring  beam  were  calibrated  as  outlined  in  Appendix  A.  Data 
were  collected  for  these  displacements  using  the  following 
procedure : 

1.  The  elastomer  was  placed  on  the  blade  root. 

2.  The  blade  was  then  aligned  inside  the  blade  mount 
and  raised  into  place  using  the  set  screw  on  the  spring  beam. 
The  set  screw  was  adjusted  to  give  the  proper  thickness  as 
monitored  on  the  microscope.  When  this  thickness  was  obtained, 
the  strain  gauge  voltage  measurements  were  recorded  and  from 
this,  the  contacting  force  was  determined  and  the  contacting 
pressure  was  calculated.  The  accelerometer  and  the  shaker 
linkage  were  then  connected  to  the  blade  tip. 

3.  The  temperature  at  the  blade  root  was  recorded  using 
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e  thermocouple  and  the  temperature  display. 

1.  The  first  displacement  fcr  a  data  run  w as 
at  the  no  slip  condition  at  t: 
was  upheld,  a  relatively  small  m.a: 
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2-t  x  cr  =  Accel- 
6.  Using  the  oscillator  amplitude 
the  desired  acceleration  was  set 
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the  blade  root  and 


liferent  blade  model  was  performed, 
r.icue  were  developed  in  a  previous 
ha  el  Kimberli.og  'Ref  13)-  This  modified 
il  root  of  one  inch  diameter  and  it  did 
inch  clearances  on  each  side.  Essentially 
viscoelastic  material  sandwiched  between 
the  mount,  the  root  fit  snugly  into  the 
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mount  and  the  damping  present  was  due  to  friction.  A  fre¬ 
quency  sweep  between  cps  and  3cC  eps  was  performed  in 

increments  of  fen  cps.  The  amplitude  of  vibrations  was  set 
to  ICCp  inches  and  the  clamp i.ng  force  was  five  pounds. 


Tata  and  Discussion 

The  data  that  were  collected  are  presented  as  plots  of 
frequency  of  vibrations  versus  force  required  to  keep  the 
amplitude  of  vibrations  constant.  The  frequency  versus  force 
divided  by  frequency  is  also  plotted  for  reasons  given  below. 

In  theory,  as  a  system  approaches  vibrations  near  a  resonant 
frequency,  the  force  required  to  maintain  a  constant  amplitude 
approaches  zero.  In  actuality,  of  course,  all  physical  systems 
possess  some  damping;  therefore,  the  force  only  drops  to  some 
minimum  value  at  resonance. 

Two  important  phenomena  are  important  to  discern  in  these 
plots.  The  first,  as  just  mentioned,  is  the  decrease  in  force 
at  resonant  frequencies.  On  the  plots  of  frequency  versus 
force,  these  points  are  easily  recognized. 

The  second  phenomenon  to  observe  is  the  flatness  of  the 
response  over  a  given  bandwidth  near  a  resonant  frequency. 


■»  vs'  Ua5  T  ^"5  ^  VO  *“  v*  ^ 
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probable  resonant  frequencies  in  a  relatively  large  range  of 
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elastomer  samples.  Analytical  predictions  of  the  resonant 
frequencies  are  discussed  in  section  III.  Values  for  the 
first  three  resonant  frequencies  were  found  to  be  26.23  ops, 
258.3  ops,  and  I"  6  4 . 9  cps  respectively.  The  first  analytically 
predicted  resonance  was  not  discernible  even  with  smaller 
frequency  increments  of  3  cps  in  this  range.  For  both  data 
runs,  there  was  an  apparent  second  resonant  frequency  at 
approximately  255-260  cps.  The  third  resonant  frequency  was 
measured  at  ’’75  ops  and  350  cps.  The  reproducibility  of  the 
data  for  each  data  run  was  quite  good  for  the  second  resonant 
frequency  and  only  fair  for  the  third  resonant  frequency. 
Figure  10  summarizes  these  broadband  sweeps  plotted  on  a 
lcgrithmic  frequency  scale  and  a  logrithmic  scale  of  Force 
divided  by  frequency  times  displacement. 

The  region  near  the  second  resonant  frequency  of  233-260 
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cps  was  chosen  to  investigate  in  greater  detail  by  performing 
frequency  sweeps  with  increments  of  5  cps. 


Figures  11,  13,  15  and  17  display  the  characteristic 
dip  in  force  at  resonant  frequencies  of  255-260  cps  for  dis¬ 
placements  of  51a  inches  and  150u  inches.  Furthermore, 
figures  12,  14,  IS  and  1?  show  the  relative  flatness  of  these 
responses  at  resonances  for  the  same  displacements. 

In  addition,  figures  19  and  2C  are  plots  of  responses 
using  friction  damping  obtained  at  frequency  increments  of 
ten  cps.  These  plots  are  presented  in  order  to  make  a  quali¬ 
tative  comparison  of  viscoelastic  damping  to  friction  damping. 

*'*  4  ’  '  -v>  Ck  ”,  **  •“  ><i  /*>  >  i  V',  "1  *"'  •*  •  -v>  Ck  T  PC  v>  y 
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the  visccelastically  damned  responses  are  shewn  to  be  scm.e- 
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sonant  irequency  of  tne  viscoelastic  system  was  arc  cps.  ...i 
higher  resonant  frequency  for  the  friction  case  is  net  sur¬ 
prising  since  this  is  a  stiffen  system  than  the  one  with  the 
elastomer  sandwiched  in  at  the  blade  root. 

The  following  is  a  list  of  other  pertinent  data  that 
were  collected  during  each  data  run: 

Run  1  Run  2 

Average  elastomer  thickness  .00414  inches  .0050  inches 

Temperature  of  blade  root  77°F  76°F 

Contact  force  at  blade  root  7.2  lbs  4  lbs 

Contact  pressure  at  blade  root  1.8  Ibs/in^  1  lb/in2 

The  area  used  to  calculate  the  above  contact  pressure  was 

determined  using  the  arc  length  of  the  dovetail  minus  the 
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BROADBAND  FREQUENCY  SWEEP 
(FREQUENCY  VERSUS  FORCE  RUN  1) 
FIGURE  6 
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BROADBAND  FREQUENCY  SWEEP 
(FREQUENCY'  VERSUS  PORCE/FREQUENCY  RUN1) 
FIGURE  7 
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Increasing  Frequency 


STIFFNESS  AND  DAMPING  '/ERSUS  FREQUENCY 
FIGURE  21 


the  fractional  derivative 


blade,  '.see  Fig  -},  the  moment  at  that  point  must  be  examined. 
Since  the  damping  material  will  be  stressed  with  a  shearing 
force,  the  above  tensile  modulus  must  be  changed  to  the 
shearing  modulus.  Then,  as  a  result. 


_  _  q  3  a  -  1_A  li 

t  t  3X 


where 


(3) 


-■2-  =  alone  of  the  beam  at 

3  x 


the  root  junction,  and 


(4) 


3  =  31C1  +  in) 

2 

=  complex  shear  modulus  (lb /in  ) 
consequently  the  moment  is  r 

M  =  t(R  gi  W)H=  a  W  t 

=  |  31(1  +  in)  ||-  (See  Fig  4) 

where 

R  =  radius  of  blade  root  (inches) 
c  =  2R-  a  V  ( inches  ) 

W  =  width  of  blade  (inches) 


(5) 


(6) 

(7) 


(3) 


Beam  Sanation 

The  classical  cantilever  beam  equation  is  as  follows 
(Ref  3) : 


(9) 


i 


■xc. 


3^  cosh  3x)  sin  3  at 


whcr51? 


jj  =  frequency  of  the  response  in  radians  per  second,  and 
8“  *  j/a  (11) 


Boundary  Conditions 

Figure  22  delineates  the  boundary  conditions  of  a  modi¬ 


fied  cantilever  beam.  The  bottom  mass  and  stiffness  at  the 


tip  represent  the  shaker  and  part  of  the  linkage  and  force 
gauge.  This  is  the  source  of  the  sinusoidal  displacement. 
The  top  mass  represents  the  remainder  of  the  force  gauge, 
linkage  and  accelerometer.  The  torsional  spring  and  viscous 
damper  at  the  root  model  the  moment  produced  by  the  shearing 
force  of  the  elastomer. 
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Consequently,  the  boundary  conditions  are: 


From  boundary  condition  (1) 


From  boundary  condition  [2) 

EI(-B2S2 3  +  3432  )  =  h  3/1  +  i'  |  3/  +  B 

let  ?  -  |  3/1  +  in) 

and  substituting  B^  =  -B2,  then 


?  2  2 
EIC-B/  sinSL  -  B/  cosSL  +  3^6  sinhsL  + 

2 

343  coshSL,;  =  3 


Using  previous  substitutions  and  rearranging 

B-,  rrL  (cosSL  +  coshSL )  +  sinhSL 

3  2EI 3 


p 

singL  -  2^r~S  (°os3-  +  coshBL) 

let  3  =  rr— v—  (cosBL  +  coshBL)  +  sinhBL 
ft 

P 

let  R  =  sinBL  -  CcossL  +  coshBL) 

'El  B 

then  3^  =  B1(R/Q) 

From  boundary  condition  (4) 

2l|^  (L,t )  -  Ma  y(L,t)  -  ks  y(L,t)  =  0 

••  ?  2 

where  yr  =  -yT ( 3  a)  ,  then 

EIa-3^30  cos8L  +  B23  sinBL  +  3,3  coshBL  + 

3 B  ^  sinheL)  +  M  (32a^2  -  k  "|  (E-.  sinBL  + 

4  Is'  S  J  1 

3 2  cosSL  +  B,  sinhBL  +;3^  coshBL)  =  0 

Using  previous  substitutions  yields 

EK-3^3'5  cosBL  -  gglg-  (B-,  +  3^(R/Q))83  sinBL  + 

31CR/Q)63  coshBL  +  (Bj  +  31(R/Q))B3  sinhBL) 

+  Ms  :s2a)2  -  ks  (B1  sinBL  -  (31  +  3,(R/Q)) 

cosBL  +  B^R/Q)  sinhBL  +  (3,  +  3  (R/Q)) 

coshBL)  *  0 

40 


V>  ■ 

*  -  O  jL 


+ 


1 

+  R/0  '  cosh 31- 1 

J 


This  problem  was  solved  using  an  iterative  technique 
delineated  in  Appendix  3.  Essentially,  the  technique  consists 
of  estimating  a  range  of  values  for  the  real  part  of  the 
eigenvalue  and  then  incrementing  the  real  part  throughout  this 
range  searching  for  values  of  the  function  that  are  decreasing. 
When  the  function  is  minimized  (ideaiy  it  would  be  equal  to 
zero)  the  resulting  real  part  of  the  eigenvalue  is  used  to 
repeat  the  above  procedure  searching  for  the  real  and  ima¬ 
ginary  parts  of  the  roots  to  the  equation.  In  this  case, 
the  search  range  is  determined  for  the  complex  part  of  the 
root . 

One  difficulty  with  this  approach  is  that  the  resulting 
minimum  value  over  a  given  range  may  not  be  an  actual  root 
of  the  equation.  To  overcome  this  difficulty,  an  order  of 
magnitude  analysis  of  each  minimum,  along  with  a  comparison 
of  eigenvalues  of  a  pure  cantilever  beam,  can  be  used  to 
eliminate  false  roots. 


There  ire  earned  programs  available  that  will  iterate  roots 
of  an  equation.  However,  since  the  roots  were  complex,  i.e.  two 
dimensional,  these  programs  were  not  used. 

Results 

The  following  were  the  results  of  the  above  program  for  the 
c  3  '_i  zi  3.  r*  *  *  cr^s  z. y-1  ^  o  o  i  • 

•  L2  y  .  Z  Z  ~  2  s  lr*j  7  c  n  c  1 0  x  F  0  o  .  x2o.2  2.-)ccSj3-.:s^^Dsi)n:=,5 


■//m; -omp-ex  Freq*  ^ ; 


ps  ;  j-,  *psi ;  n=  .c 

-  - j.uS-  ,ri  —  •  -  — 


3  -  *  (  .  0  ?  9  -  ,  .  11”  )  /  in ;  lamp  lex  Free.  (‘'64.9>37.4)ccs;3,*2l4psi;  n  =  l.l’7 

In  this  case,  the  mass  and  stiffness  of  the  shaker  were  removed 

along  with  part  of  the  linkage  mass.  The  remaining  mass  was 

0.-0ClJ3  los  -  sec  /in  ,  attributed  to  the  linkage  past  the  force 

gauge,  one  half  the  force  gauge  mass  and  the  accelerometer  mass. 

The  reason  for  including  this  mass  in  the  configuration  is  that 

at  resonance,  the  force  should  be  at  zero.  Consequently,  the 

mass  and  stiffness  of  the  shaker  along  with  part  of  the  linkage 

was  removed  from  the  circuit.  With  these  eigenvalues,  the  re- 

0 

suiting  frequencies  were  calculated  using  equation  11  and  a/1  = 

w .  The  following  is  a  list  of  the  values  of  the  other  parameters 

that  entered  into  the  equation  for  the  frequencies: 

=  ( .5363/.3C4l4)(U3  +  i  21.5)  (10-5)  =  (1.-,6  +  t  ,  for 
21  31(12)°' .001373)  '  1 

PL  =  1-  6363/,  0^l4)(133  +.1  lr-)(i:o)  ,  :5.55  +  i  5.2o)  for  8? 

31(10)  (.101373) 

?L  .  (-5369/.  3C.ll4)(2l4.t.l. 250)  .  (3>?5  +  ,  10-2U)  for  s 

ij-  31(10)'J( .  001373) 


Ms  a  ~  _  Ms 

Zj  j.  u  Mfc  - 


1 .  0001- 


1 4  0  4  }  (  1 0  . 


0 


k.  L3 
s 


0  C1Q.5I3 _ 

31(10)°  (0.001373). 


A  comparison  of  the  analytical  and  experimental  results  is  as 
follows : 


Resonant 

Analytical 

Experimental 

Frequency 

Prediction 

Measurement 

3  -  **  vi  0  rj  -i  ^ 

Second  Run  1 

253 

255 

1 . 1 

Second  Run  2 

253 

260 

negligible 

Third  Run  1 

764 

775 

1.4 

Third  Run  2 

714 

353 

11 

.T f.  i  3o  iy  Mo 

of  the  lowest  mo 

de ,  e mrected  to  be 

d  r*0  d  cr>  insnt 

'‘  nm'  t?  t*  —  -  3  r  v~  vn';  0  rr 

vi  a  t;  n  vp  2“ 

Ut  '  O'" 


Oyi  *"  ^  v»  - 


—  i  n )  -? 


is  the  mass  moment  of  inertia  about  the  blade  dovetail 
root  made  up  of  the  mass  of  the  blade  and  the  tip  mass  des¬ 
cribed  in  the  previous  section.  This  tip  mass  is  treated  as 
a  point  mass;  hence,  the  mass  moment  of  inertia  about  its  own 

center  of  rravitv  is  ignored.  In  short,  I  is  calculated  as 
^  J  o 

follows : 


i+3 


-  J  V3  +  VLl2 

=  -  C0.0Q46)(10.5)3  +  0.000143C10.5)2 
=  0.1348  lb  -  sec2  -  in 

Here  the  moment-  of  inertia  of  the  dovetail  has  been  ignored. 

The  moment  at  the  blade  root  is  similar  to  equation  6 
but  for  the  rigid  body  mode  ■—£.  j_3  replaced  with  •> .  Hence, 
the  real  part  of  the  natural  frequency  for  equation  (.12).  is 
calculated  as  follows: 

-»  n 

2  w 

jj  =  — — — 

***  3  " 

3-*  =  vi ^  b  / ^  n**  n  "3  •*■*•'»  ^  <~*  *  *  o ^ ^  ^ o. ^ ^ ^  *- v» ^ 

frequency  is  as  follows: 

_  ’(.5363)03)  1  l 

^  ('  “  Y7  ”  ^  “  ’I  t 

=  136.5  rad/sec 
=  31.25  cps 

Since  this  frequency  is  close  to  the  one  predicted  in  the  beam 
equation  analysis  -26.2  cps  -  it  does  not  explain  why  it  could 
not  be  found  experimentally  as  a  resonance  . 

As  a  final  attempt  to  explain  this  enigma,  an  analytical 
prediction  of  the  response  was  examined.  The  hypothesis  was 
that  if  the  damping  in  the  rigid  body  mode  was  sufficiently 
heavy,  then  the  dip  in  force  at  resonance  may  not  be  discern¬ 
ible.  The  analysis  showed  that  this  resonance  should  de- 
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period  reouired  to  ;ol2ect  one  data  run,  one  thermcccucle 


reading  was 

20  gP'”  1?  OQ  “  VQf’V  O 

oo  i:3?. 

T?  ■»  7’  1  Vlp  g 
~ 

,  C  v*»  ri  ^  3  V>i  f  *“  *■*  C. 

r*«  d"*  *-  ->,r» 

V.  «  fl  3  ^ 

tioned  above 

•”nV\^v»0  o 

S  r  ^*''ov»or'r 

.oss  factors  clotted 

to  account  f 

or  the  temperature 

variat ions 

for  both  the  second 

and  third  resonant  frequencies. 

The  data 

used  for  this  ccmpar 

ison  are  as 

follows : 

Loss 

Modulus 

Temperature  Frequency 

Factor 

( lbs/in2 ) 

°C 

(CPS) 

n]_  =  1.1203 

138 

24.4 

315.1 

n2  =  1.0131 

117 

29-4 

304 . 4 

=  .8981 

78 

37.3 

284.3 

In  addition,  the  elastomer  thickness  was  varied  from  0.003 
to  0.005  inches.  The  abcissa  is  scaled  to  the  magnitude  of 
the  complex  modulus  divided  by  the  thickness. 

Figure  25  shows  that  over  the  entire  range  of  tempera¬ 
tures  and  thicknesses  the  predicted  resulting  frequency  vari¬ 
ation  is  less  than  13%.  Figure  26  shows  that  at  the  higher 
frequencies  this  error  is  reduced  to  less  than  7%.  In  sum¬ 
mary,  variations  in  the  above  parameters  will  probably  not 
cause  large  discrepancies  in  the  results  of  the  experiment. 
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V  Conclusion 


In  conclusion,  viscoelastic  materials  appear  promising 
as  an  effective  vibration  damper  for  turbine  and  compressor 
blades.  Presently,  most  turbine  engines  use  friction  damping 
at  the  blade  root  of  the  turbine  and  compressor  blades.  3y 
qualitatively  comparing  the  flatness  of  response  curves  at 
resonance  for  friction  and  viscoelastic  damping,  the  latter 
was  found  to  have  a  flatter  response.  Hence,  it  appears  to 
be  a  more  effective  damping  mechanism.  However,  before  this 
application  can  come  to  fruition,  more  research  must  be  per¬ 
formed.  Important  questions  must  be  answered.  Questions  such 
as:  (1)  "Is  the  increased  complexity  of  a  viscoelastic  blade 

roo1:  material  overcome  by  the  benefits  of  improved  vibration 
damping?''  and  (2)  "Can  a  design  be  formulated  which  is  ser- 
vicable  in  the  adverse  environment?" 

Finally,  this  was  intended  zo  be  a  preliminary  study. 
During  the  course  of  the  study  three  other  observations  were 
made  as  follows: 

(.1)  holding  displacement  constant  while  measuring  the 
force  is  a  feasible  approach 

(2)  even  though  the  analytical  model  is  sensitive  to  tip 
mass,  such  a  model  can  be  made  that  correlates  with  experi¬ 
mental  data  at  selected  points 

(.3)  if  higher  displacements  had  been  used  at  the  first 
resonant  frequency  the  force  measurements  may  have  been  greate 
than  the  noise  level. 
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STRAIN  GAUGE  CALIBRATION  CURVE 
FIGURE  27 


V  1.3  3  3  4-  ■**'  ■£***  &  g  s 

F  is  a  known  force 
D  is  a  displacement 

The  procedure  consisted  of  positioning  the  shaker  so 
that  the  core  motion  was  vertical.  A  known  dead  weight  was 
balanced  on  the  tip  of  the  element  and  the  displacement  was 
measured  with  the  Gaertner  microscope,  so  that  a  stiffness 
could  be  computed. 

Next,  an  accelerometer  was  mounted  on  the  tip  of  the 
element.  Using  an  appropriate  amplifier,  the  accelerometer 
output  was  connected  to  a  brush  recorder  model  RD  2522  20. 
The  shaker  was  given  a  brisk  impulse  and  the  output  of  the 
oscillations  were  recorded  on  the  Brush  recorder.  With  the 
recorded  calibrated  time  axis,  the  fundamental  frequency  of 
oscillation  was  determined.  Then  using  using  the  formula: 


5° 


where 


oj  is  the  natural  frequency 

k  is  the  stiffness  determined  above 

M  '  is  the  mass  of  the  shaker,  connecting 
bolt  and  accelerometer 

the  moving  mass  of  the  shaker  with  the  accelerometer  and 
connecting  bolt  was  calculated  from  the  formula  and  the  mea¬ 
sured  stiffness.  The  accelerometer  was  weighed  and  its  mass 
was  substracted  from  the  calculated  total.  Since  the  con¬ 
necting  bolt  remained  on  the  shaker  for  all  the  experimental 
tests,  its  mass  was  not  subtracted  from  the  total.  The 
resulting  values  are  as  follows: 

k  =  11.54  lb/in 

u)n  =  11.5  cps  =  72.26  rad/sec 

M  =  M  '  -  accelerometer  mass 
s  s 

=  0.002048  lb-sec^/in 


any. 

In  the  first  half  of  the  program,  the  imaginary  parts 
of  the  complex  arguments  are  set  to  zero  to  obtain  a  rough 
estimate  of  the  real  part  of  the  eigenvalue.  With  this  es¬ 
timate,  the  second  part  of  the  program  iterates  on  both  the 
real  and  imaginary  parts  of  the  actual  eigenvalue.  Then,  the 
range  of  the  eigenvalue  estimate  is  increased  and  the  above 
procedure  is  repeated  for  the  next  eigenvalue. 


PROGRAM  EQROOT 

C  THIS  PROGRAM  IS  DESIGNED  TO  ITERATE  TO  SOLUTIONS  OF  THE 
C  EIGENVALUE  PROBLEM  RESULTING  FROM  THE  EQUATION  OF 
C  MOTION  IN  THE  ANALYTICAL  SECTION 

REAL  G1,G2,C,T,L,M,K,J,F1,F2,F3,F4,X1,X2 
COMPLEX  Z,Z1,B,D,F,W 
COMMON  G1,G2,C,T,L,M,K,J,A,Z1,W 
Z1«(0. 0,1.0) 

C»(2. 5)* (0. 496**3.0)* (2. 2514) 

C  C  IS  THE  CONSTANT  FROM  EQUATION  7 
L*10 . 5 

C  L  IS  THE  LENGTH  OF  THE  BLADE 
M«0. 002458 

C  M  IS  THE  MASS  OF  THE  SHAKER,  LINKAGE,  FORCE  GAUGE  AND 
C  ACCELEROMETER 

K-11.5^ 

C  K  IS  THE  STIFFNESS  OF  THE  SHAKER 

J-31*( 10. 0**6)* (0.001373) 
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-  n  myr  op/^nr/iw 


:  5  IS  THE  COMPLEX  LOWER  RANGE  0?  THE  EIGENVALUE 

RE AO*,  3 
PRINT*  ' XI®' 

G  XI  IS  THE  UPPER  LIMIT  OF  THE  REAL  PART  OF  THE  EIOENVALU 
READ*, XI 

13  I ?( RE AL ( 3 ). LE . X 1 } THEN 
3=3+D 

?2=REAL(F(3 }  } 

IF(F1**2.0.GE.?2**2) THEN 
GO  TO  10 
END  IF 
3*B-2*D 

IF  (REAL(D) .LT. 0.0001) THEN 
GO  TO  20 
END  IF 
D-D/ ( 10 ) 

GO  TO  10 
20  PRINT*,’ 3=’ 

PRINT*, 3 

W-((A»3**2.0)/(2.0*3.1U16)) 

PRINT*, »W*’ ,W 
PRINT*, 1 Fl= ’ ,?1 
PRINT*, *P2«» ,F2 
END  IF 

0  THE  ESTIMATE  OF  THE  REAL  PART  OF  THE  EIGENVALUE 

C  IS  USED  TO  ITERATE  3EL0W  ON  THE  COMPLEX  ARGUMENT  OF 

C  THE  EIGENVALUE 

PRINT*, ’3=?’ 

RE AD*, 3 
PRINT*, ’D-» 

READ* ,D 
PRINT* , ’ X2* ’ 

READ* ,X2 

C  X2  IS  THE  UPPER  RANGE  OF  THE  COMPLEX  PART  OF  THE 

0  EIGENVALUE 

30  IF ( AIMAG ( 3 ) . LE . X2 ) THEN 
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3=3+?.EAl 3 

f£=aimag-;f'3) } 

IF- ?1**2.  I+F2**2.0.  SE.F;**2.  0+Fl**2.  :  the: I 
30  TO  -10 
END  IF 

3=B-2 . 0* ( REAL(D ) ) 

D-D/10 

-  V  j'J 

50  PRINT*  ' 3= f  3 

W-(A*B**2. 0)7(2. 0*3. 1416) 

PRINT* , 'W* '  ,W 
PRINT*, ' FI* ' ,F1 
PRINT* , ’  F2-’ ,F2 
PRINT*, 'F3=' ,F3 
PRINT* , '  F4-  ’ ,F4 
END  IF 
END 

COMPLEX  FUNCTION  ?(Z) 

THIS  FUNCTION  COMPUTES  A  VALUE  FOR  THE  EQUATION  OF 

MOTION  THAT  SHOULD  CONVERGE  TO  ZERO  TO  SOLVE  THE 

EIGENVALUE  PROBLEM 

COMMON  G1,G2,C,T,L,M,X, J,A,Z1,W 

COMPLEX  Z,Z1,3,W,Y 

REAL  G1,G2,C,T,L,M,X,J,A,D,X 

COMPLEX  P ,R,Q 

Z  =  Z*L 

?=(C/T)*(G1+Z1*G2) 

R=CSIN( Z ) - ( CCOS ( Z ) +CCOSH ( Z ) )*(?/( 2 . 0*J* ( Z/L) ) ) 
Q=CSINH(Z)  +  (CCOS(Z)+CCOSH(Z) ) * ( P/( 2 . 0  * J* ( Z/L) ) ) 
F=-CCOS ( Z ) - (P/( 2 . 0  * J* ( Z/L ) ) ) * ( 1 . 0+R/Q ) *CSIN ( Z ) 

*+(  R/Q)  *CCOSH(Z  )  +  (?/(_2.0*J*(Z/L)  )  )*(  1 . 0+R/Q)  *CSINH(  Z ) 
»+(1.0/CJ*(Z/L)**3.0) )*(M*(A*(Z/L)**2.0)**2.O 
*-K) *(CSIN(Z)-(P/(2.0*J*(Z/L) ) )*(1. 0+R/Q) *CCOS(Z)+ 
*(R/Q) *CSINH(Z ) +( P/(2 . 0*J*( Z/L) ))*( 1 . 0+R/Q) *CCOSH( Z ) ) 
Z  =  Z/L 
END 

COMPLEX  FUNCTION  CSINH(Z) 

THIS  FUNCTION  COMPUTES  A  VALUE  FOR  THE  COMPLEX 
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Response  Analysis  Rigid  R° rZ 

The  equation  of  motion  for  a  rigid  body  with  viscous 
damping  and  harmonic  excitation  is  as  m-lows  .Ref 

•nvCt'  +  kd  +  1y)  y(^'  =  AkelU>t 


Assuming 


mm 


iujt  _  g  COs  tut  +  Bi  sin  it 


-oT  Bmeiu)t  +  k(l  +  it)  Selu)t  =  Ake"w 


B  *  - - - T~ 

k(l  +  iy)  m 


1  +  iy  (m/k) 


Hi  *  i— _  }  therefore 


l  +  iy  -  — r 

V 


The  amp 


litude  of  the  steady  state  displacement  is 


then 


.  .i  c  w  \  2  i 

1  +  iy  ; 

n 


.  .  >2,2  2n1/2 
(«/»n)  )  +  Y  ' 


1 

1 


In  order  to  match  the  physics  of  the  turbine  blade  pro¬ 
blem,  equation  (13)  should  be 

I  $  +  3,  ( c/t ) ( 1  +  iy)  *  *  M(t)  =  P  Lelut 

O  j_  o 

but,  the  following  parameters  can  be  substituted: 


y  -*■  <pL 


L  (c/t) 


m 


Ak  =  F 


Y  = 


:cnsequently ,  equation  (14)  can  be  written  as 


31(c/t) 


, .  2,  2.2  ,  2 

( 1  -  w  /u>  )  +  n 


1/2 


|y 


(15) 


Figure  30  is  a  plot  of  FQ  versus  w  for  frequencies  be¬ 
tween  10  cps  and  50  cps.  The  following  values  were  used  in 
equation  15  when  determining  the  values  of  Fq : 

Gx  =  43  lb/in2 

x  =  150  u  inches 

n  =  .5 

U)  =  31.26  cps 


t 


r 


> 


a 
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Analytical  Prediction  of  Force 
Versus  Frequency  for  Rigid  Body 
Mode  at  150  / JL  Inch  Displacement 
Figure  30 
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